Chiral symmetry breaking and 
quark confinement 
in the nilpotency expansion of QCD 

Sergio Caracciolo 
Dipartimento di Fisica and INFN 
Universitd degli Studi di Milano 
via Celoria 16, 1-20133 Milano, ITALY 
Sergio . CaraccioloOmi . infn. it, 

Fabrizio Palumbo 
INFN - Lahoratori Nazionali di Frascati 
P. O. Box 13, 1-00044 Frascati, ITALY 
fabrizio .palumboOlnf . infn. it 



June 6, 2011 

Abstract 

We apply to lattice QCD a bosonization method previously devel- 
oped in which dynamical bosons are generated by time-dependent Bo- 
goliubov transformations. The transformed action can be studied by 
an expansion in the inverse of the nilpotency index, which is the num- 
ber of fermionic states in the structure function of composite bosons. 
When this number diverges the model is solved by the saddle-point 
method which has a variational interpretation. We give a stationary 
covariant solution for a background matter field whose fluctuations 
describe mesons. In the saddle-point approximation fermionic quasi- 
particles exist which have quark quantum numbers. They are confined 
in the sense that they propagate only in pointlike color singlets. Con- 
ditions for chiral symmetry breaking are determined, to be studied 
numerically, and a derivation of a mesons-nucleons action is outlined. 



1 Introduction 



The fundamental fields appearing in QCD, that is quarks and gluons, are 
confined, and we are able to observe directly only mesons and baryons. It is 
conceptually interesting to see how these composite fields emerge from the 
microscopic dynamics and in several cases can be practically convenient to 
reformulate QCD in terms of them. The description of the phase diagram 
of the theory, in particular, should be more transparent in terms of these 
fields. We will refer to a reformulation of QCD in terms of hadronic fields as 
to QCD hadronization. 

This task goes beyond the perturbation theory in the gauge coupling 
constant, thus our starting point will be the lattice formulation of QCD yy. 

After Wick rotation, in the Euclidean formulation, the path-integral of 
the pure-gauge sector could be studied using a number of tools, first of all 
Monte Carlo simulations, which have helped us to understand what occurs 
in the non-perturbative regime. Within this framework, the role of quark 
fields soon appeared to be more difficult to consider. There is a funda- 
mental difficulty with a lattice action for fermions which explicitly preserves 
locality, gauge and chiral symmetries: the duplication of the spectrum (see 
for example the textbook [2]). But also, as in the path-integral formulation 
fermionic fields are represented by Grassmann variables in a Berezin integral, 
an efficient numerical simulation seems to require the preliminary integration 
of the fermionic degrees of freedom. The recovered functional determinant 
is heavily demanding from the numerical point of view and is still a hard 
problem in the region of finite chemical potential because it is not positive 
definite. This prevents the definition of a probability measure, thus making 
the introduction of approximation schemes difficult. 

A hadronization of QCD could overcome these difficulties, at least in the 
mesonic sector. We have, therefore, been pushed to apply to lattice QCD a 
general method of bosonization we developed in recent years, both in Many 
Body and Field Theory, in the presence of fundamental fermions. It turns 
out that not only does such a method allow us to introduce mesons, but also, 
in some approximations, baryons. 

Our starting point is the operator formulation of the partition function in 
the Fock space representation for the fermionic fields, in which approxima- 
tions can be introduced by following physical insight and can be mathemat- 
ically justified. In a first approach [2111] we restricted the evaluation of the 
partition function of a system to states of composite bosons, in a variational 



2 



spirit. The resulting bosonic effective action in such restricted space was 
evaluated exactly. 

In such a framework a perturbative scheme was formulated as follows. 
The composite bosons are characterized by an integer, the index of nilpo- 
tency, which is the number of fermionic states in their structure function. 
This is also the maximum number of composites which can exist in a given 
quantum state and therefore only when this number diverges can composite 
bosons behave as canonical bosons [Ij. We set a perturbative expansion in 
the inverse of this number which we call nilpotency expansion. Since the 
index of nilpotency counts the total number of independent fermionic modes 
in the composites, it is in general much greater than the number of inter- 
nal degrees of freedom of the fermions. For instance the number of degrees 
of freedom of the electron is 2, but the total number of fermionic states in 
the Cooper pairs in superconductors is infinite in the thermodynamic limit, 
which is the reason why the BCS solution gives the exact energy per particle 
in this limit. 

It was later recognized [5] that this variational approach provides the 
saddle-point approximation in the nilpotency expansion to a theory obtained 
from the original one by time-dependent Bogoliubov transformations. Bogoli- 
ubov transformations are a natural way to introduce composites as Cooper 
pairs (see for instance [6] for a different approach along similar ideas, and [7] 
for the solution, directly in the continuum, of the problem of Dirac parti- 
cles in an external stationary magnetic field in 2+1 dimension by means of 
a series of Bogoliubov transformations). But the specific difficulties of the 
renormalization procedure in the Hamiltonian formalism have limited their 
use in this domain. A more severe drawback in an application to gauge the- 
ories is that symmetric terms of the original theory, in general, give rise to a 
sum of terms in the transformed actions none of which conserves in general 
the symmetries (even though, since Bogoliubov transformations are unitary, 
all symmetries are conserved in the sum). This becomes a potential source 
of problems when some approximations are performed. 

In our approach we try to avoid difficulties with renormalization in the 
Hamiltonian formalism by use of the transfer matrix formalism, which is 
at our disposal since we are using the lattice formulation. By using this 
formalism we can oscillate between the operator formalism in Fock space 
and the functional formalism. 

An independent Bogoliubov transformation at each time slice is per- 
formed in the operator form of the partition function. The time-dependent 
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parameters of the transformation are required to vary under symmetry trans- 
formations in such a way that the quasiparticles fields transform in the same 
way as the original fermion fields under the symmetries of the theory, in par- 
ticular gauge invariance. These parameters can then be associated with dy- 
namical bosonic (composite) fields in the presence of fermionic fields (quasi- 
particles) with the quantum numbers of the bare fermions. A compositeness 
condition avoids double counting [S]. One thus gets an effective action of 
composite fields plus quasiparticles, exactly equivalent to the original one, 
in which ground and excited states can be treated on the same footing. Of 
course, in practical applications, some approximation must be introduced. 

As usual different solutions in a saddle-point approximation can be re- 
lated to different phases of the theory and the nilpotency expansion will 
be our tool to study the dynamics of the composites. There is a complete 
arbitrariness regarding the composites introduced by the Bogoliubov trans- 
formations because, since they are unitary, the transformed theory is exactly 
equivalent to the original one irrespective of their choice. But only when 
the introduced composites reproduce the effective degrees of freedom in the 
given phase, the transformed theory, after suitable approximations, can have 
practical applications. 

In a successive investigation [8] of our method, we looked into the na- 
ture of the saddle-point equations. In the saddle-point approximation the 
Bogoliubov transformation on the transfer matrix has the same effect as the 
Foldy-Wouthuysen transformation on the Dirac Hamiltonian: it eliminates 
the direct mixing between fermions and antifermions. In the absence of gauge 
fields we found an explicit solution of the saddle-point equations ^4]. 

In the present work we consider the saddle-point equations in the presence 
of gauge fields, namely the application of our method to QCD, restricting 
ourselves to the case of zero chemical potential. A preliminary discussion 
of finite baryon density can be found in [9]. We use the Kogut-Susskind 
regularization for fermions in the flavour basis. The corresponding expres- 
sions in the spin-diagonal basis which is more commonly used in numerical 
simulations will be reported elsewhere. 

We find an exact covariant solution which requires that the vacuum should 
be dominated by stationary chromomagnetic fields. At the present stage of 
our research such a dominance appears to be driven by dynamical quark 
fields, but it should be remembered that dominance of chromomagnetic fields 
was already found and discussed by several authors for pure-gauge theo- 
ries [inilll]. Because of the coupling with quarks, the QCD vacuum appears 
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to be a dual superconductor (not a color superconductor). Dominance of 
cliromomagnetic fields in the non-Abelian case, is indeed reminiscent of the 
picture of color confinement based on electric-magnetic duality p^VflG] and 
has been found also in the standard theory of color superconductivity [T7] . 
This should be compared with the abelian case, where, by contrast, we do 
not expect that the configurations of gauge fields giving rise to a stationary 
nonvanishing magnetic field should dominate the vacuum. 

Our vacuum contains a condensate of quark-antiquark pairs, but these 
Cooper pairs do not have the quantum numbers of chiral fields, which are 
instead associated with the fluctuations of this condensate. For this reason 
we will refer to this condensate as to the background field. 

Does this background solution describe what we expect as essential fea- 
tures of QCD like confinement and chiral symmetry breaking? Our answer 
can be only partially positive. Concerning the first issue we obtained an 
important result: in the background of a stationary chromomagnetic field 
quasiparticles cannot propagate separately: only color singlet composites of 
quasiparticles can propagate and therefore have a physical particle interpreta- 
tion. This means that color is confined as far as quasiparticles are concerned. 
The quasiparticles color singlets can be mesons and baryons as well. There- 
fore not only can our method describe bosonization, it can also account for 
the more complicated process of formation of composites of an odd number of 
quasiparticles. Concerning spontaneous chiral symmetry breaking instead we 
we have found conditions similar to the standard ones which must ultimately 
be solved numerically. 

The expression of the background field is the starting point of our program 
of hadronization. This will require, at the present stage, the evaluation by 
numerical simulations (possible because of our lattice regularization) of some 
quantities appearing in the nilpotency expansion. An illustrative example of 
such quantities can be found in Section [TJ where a derivation of a meson- 
nucleus action is outlined, and gauge-invariant coefficients depending on link 
variables appear explicitly. Also, the expression of the partition function at 
finite chemical potential that we found [S], which avoids the sign problem, 
requires an integration over spatial gauge link variables. 

We wish to remark that our method is well suited to also study the effect 
of an intense background magnetic field on strong interactions, a problem 
considered of interest both at the level of the cosmological electroweak phase 
transition and for the heavy-ion collisions. And, indeed, numerical simula- 
tions have already been performed both in the quenched approximation, see 
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for example pE], and with dynamical fermions (see also [IS] for a detailed bib- 
liography), in order to try to understand magnetic catalysis, i.e. the increase 
of chiral symmetry breaking induced by the magnetic background field. 

The paper is organized in the following way. In Section [2] we establish 
our definitions and notations. To make the paper self-contained we report in 
Section [3] the derivation of the effective action by time-dependent Bogoliubov 
transformations. In Section|l]we report our solution of the saddle-point equa- 
tions for gauge theories. In Section [5] we discuss the physical interpretation 
of our results and their relevance to symmetries, in particular chiral symme- 
tries. In Section |6] we prove confinement of quasiparticles in the saddle-point 
approximation. In Section [7] we outline a derivation of mesons-nucleons ac- 
tion from QCD. Finally, in Section [SI we perform a summary and give an 
outlook of our method. 

2 Definitions and notations 

Consider a system of fermions interacting with external bosonic fields includ- 
ing gauge fields, regularized on a lattice. The fermionic part of the partition 
function at finite temperature T can be written 

Lo/s-l 

t=o 

Lq = T^^ is the number of links in the temporal direction, T is the fermion 
transfer matrix, Tr^ is the trace over the Fock space of fermions. The pa- 
rameter s takes the value 1 in the Wilson formulation for lattice fermions, 
but s = 2 for the Kogut-Susskind fermions which live on blocks of twice the 
size of the lattice spacing. The index t labels the blocks along the "time" 
direction. 

For Wilson fermions the expression of T was given by Liischer [20] , in the 
gauge Uo = 1, in which one has to impose the Gauss constraint in the Hilbert 
space of the system (a Fock space of fermions in which the coefficients of the 
fermionic states are polynomials of spatial link variables). One can also use 
a slightly modified form which avoids the Gauss constraint by reinstating the 
temporal links variables: 

rt,t+i := e^'-(*^*+*^')r/ V^e^^"fi+i (2.2) 
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where n is the chemical potential and h is the fermion number operator 

n := v)u — v'^v , (2.3) 

(the sum on all fermion indices is understood) with and v"^ (respectively 
u and -0), creation (respectively annihilation) operators of fermions and an- 
tifermions, obeying canonical anticommutation relations and 

ft = exp[-u^MtU-v^M^v]exp[vNtu\ (2.4) 
Vt = exp[uhnUo,tu + vHnUltv]. (2.5) 

The matrices Mf {Mf being the transposed of Mj) and Nt are functions of 
the spatial link variables at time t and possibly of other bosonic fields, such 
as the external magnetic fields considered in [T8|IT9]. 

Remarkably, Kogut-Susskind fermions, but in the so-called flavour basis 
(see for example [2J), give rise to a transfer matrix of the same form j2T]. 
Explicit expressions for Wilson and Kogut-Susskind fermions in the flavor 
basis are reported in Appendix |Al The variables Uo^t are matrices in a uni- 
tary representation of the gauge group whose elements are the link variables 
between Euclidean time t and t + 1 

(f/o,t)x,y = 5x,yf/o(t,x) (2.6) 

where boldface letters, such as x, denote spatial coordinates. 

We introduced the following notation, which we will use for any matrix 

A 

tr±A := tr (P±A) . (2.7) 

The operators P± project on the components of the fermion field which prop- 
agate forward or backward in time 

it = P+ip 

= P-^j (2.8) 

and their expressions are given in Appendix |Al The symbol "tr" denotes 
the trace over fermion-antifermion intrinsic quantum numbers and spatial 
coordinates (but not over time). Finally we will denote by Tq^'' the forward 
and backward translation operators of one block, that is s lattice spacing, in 
the "time" direction 

[T^Xm = Kt.±i ■ (2.9) 
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The expression of the full partition function is 



Z = j [dU] expi-Sa) Zf 



(2.10) 



where Sq is the pure gluon action. 



3 Time-dependent Bogoliubov transformations 

The material of this Section is taken from Ref. [8] and is reported to make 
the paper self-contained. 

Usually the trace appearing in the definition of the transfer matrix is 
evaluated using at each time slice coherent states of fermions 



|«, l3) = exp(— a-u''" — /3v'^)\0), 



(3.1) 



where the a, (3 are Grassmann fields. We will use instead states obtained by 
applying, at each time slice, an independent Bogoliubov transformation. 



T 



where 



1 + J^J^^)^^ 



i? = (l+^t^)-i ^ = 
and T is an arbitrary matrix such that 



(3.2) 
(3.3) 

(3.4) 



The circle over the R denotes the involution defined by the above equations. 
The new operators (we omit the subscript J-" to lighten the notation) 



f3 



(3.5) 
(3.6) 



satisfy canonical commutation relations for any choice of the matrix J-". We 
will let J-" depend on all the fields coupled to the fermions in such a way as 
to respect as many symmetries as we can. The vacuum of the new operators 
is 

1^) = exp(^t) |o) (3.7) 
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where 

jrt = u^j^^\ (3.8) 

is a the creation operator of a composite boson. As aheady mentioned the 
new vacuum appears as a coherent state of fermion-antifermion pairs. We also 
remark that the new vacuum defined in (13.71) is gauge-invariant and therefore 
satisfies the Gauss constraint. The transformed states can be written 



#(J^) =(det+i?5) |a,/3;jr) (3.9) 

= (det+i?5) exp(-aa^ -/3/3'f)|J-) (3.10) 
= (det+i?5) exp i^P^ -aa^ -bP^ - pj^o^ |0) (3.11) 

where det± is the determinant in the subspace where P± projects, with 

det+i? =det_i? (3.12) 
det±i?=det±^ (3.13) 



1 ° 1 

and a := R~2 a and b := (5 R~2. The exphcit definition of the operator ^ 

can be found in [8^ Appendix B], and here we correct a misprint in [SI (2.16- 

2.18)] where the normahzation factor det^ had been forgotten. 

After evaluation of the trace, the partition function becomes 

Zf =exp{-5^,(^)} / exp{-^,p(a, /3; J")} (3.14) 

where the Grassmann variables a*,a, P*, P satisfy antiperiodic boundary con- 
ditions in time. In the above equation Sme, the term independent of the 
Grassmann variables, will be interpreted as a meson action 

Lo/s— 1 Lo/s-1 

Sme{J'):=- Yl tr+ In (i?t [/o,t = - ^ ti + In {Rt Et+i,t) (3.15) 

t=0 t=0 

where 

Et+i,t := {J'N,t+i)^ e^'^^'Ul e''' J^N,t + Pt+i e"*'*+^^o^ e"*'* J't , (3.16) 
with 

TN,t-=l + N}Tt (3.17) 
and we used the fact that Uo^t are unitary. 
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The other term is the action of quasiparticles 

Lo/s-l 

t=o 

+ aliVt - Ht)at+i - A+i(Vt - Ht) A*] (3.18) 
written in terms of lattice covariant derivatives 

Vt-.^s-'l^e'^'Uo^t-Tt^) (3.19) 
^o,t--o^) (3.20) 
and the lattice Hamiltonians, respectively, for fermions and antifermions, 

Ht := s-'e^^ (Uo,t - (3-21) 
Ht := s-'e-^^ (Ul- ET\, hf^) . (3.22) 

There are, in addition, unwanted terms which mix quasiparticles with quasiantipar- 
ticles whose coefficients are 



Rt- E7U J^N,t-i e^*-if/o,t-ie 



(3.23) 



Rt-e * Uo,te [TN,t+i) -Bt+i,; 



o 1 



Rl ■ (3.24) 



The definitions of the other new symbols are 



Et+i,t:^ h,t e^*Vo,te*^'+^ (^iv,t+i) V T"* e-^^*Vo,te-^*+i J-.V^ (3.25) 



(3.26) 



4 Saddle-point equations and factorization of 
the transfer matrix 

We assume that the contribution of quasiparticles to the vacuum energy be 
negligible. Therefore in order to determine the contribution of the fermions 
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to the vacuum energy we must minimize the mesonic action with respect to 
J-", J-'^ This gives the saddle-point equations, vahd for < t < — 1 

= iV,+i + e-^^'+^f/oV-^-^*V,(j-^,0"'e-^-^'Vo,e-^^'+^ (4.1) 
A = Nl + e-^-^'Vo,e"*^'- J-^ie-^^-^f/oV"*^'^ . (4.2) 

The main difficulty of the saddle-point equations stems from their depen- 
dence on time. This difficulty is reduced if we look for stationary solutions, 
as appropriate to the vacuum. If J-" is stationary, the elementary bosonic 
fields coupled to the fermions which enter its expression should also be sta- 
tionary [3]. In gauge theories J-" certainly depends on spatial link variables. 
Stationarity for gauge fields can be formulated in a gauge covariant way by 
requiring that these fields evolve according to gauge transformations, so we 
must require that 

f/fc(t,x) = iy,|,^,(0,x)iy,_,^fc. (4.3) 

As a consequence the chromomagnetic contribution to the pure gauge-field 
action, namely the trace of spatial plaquettes, does not depend on time. 

Accordingly, the matrices Nt, Mt are related to those at time t = 0, that 
is, if iVo = iV and Mq = M, by 

Nt = W} NWt, Mt = W} M Wt . (4.4) 

We still wish to set the contribution of the chromoelectric field to the 
gauge action, namely the trace of spatio-temporal plaquettes, to be indepen- 
dent on time. We have been able to arrive at a stationary solution for J-" 
only with the particular choice 

Wt+i,^ = f/o(0, x) f/o(l, x) . . . f/o(t, x) (4.5) 

which lets the contribution from the chromoelectric field vanish at all times. 

Indeed, if we do not consider colored composites, the saddle-point equations 
for 

Ft = W} J-o Wt = W} T Wt (4.6) 

then become 

J^ = N + e-^^e-^V J-^^e-^-^'e-*^ (4.7) 
and the Hermitian conjugate relation. 
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Assuming M = M\ a condition satisfied in both the Wilson and Kogut- 
Susskind regularizations, under the condition [A^, M] = (which is not sat- 
isfied by the Wilson regularization in presence of a nontrivial gauge configu- 
ration), their solution [UIH] is 



= N{2N^N)-^ -Y + ^Y"^ + ANm . (4.8) 

where 

Y = 1- N^N - e-^^ . (4.9) 
The time evolution of the quasiparticle Hamiltonians is slightly different 

= w} H Wt+, , nt= wl, n Wt . (4.10) 

At the saddle point (we will overline all quantities evaluated at the saddle 
point) 



1 - :F^'e-2^ 



(4.11) 



so that "H and "H are Hermitian functions of M and A^^A^ and the vacuum 
energy is 

Sme = S^T) = -— tr+lng (4.12) 
s 

where we introduced 

Q = (1 -se-^^'W)"^ (4.13) 

for future convenience. 

For Wilson fermions we have not been able to find an exact solution, but 
in the following we will assume its existence. 

Since the matrices A^, M are Hermitian and by assumption commute with 
each other, they can be diagonalized simultaneously and J-" is diagonal in such 
a basis. Labeling each eigensubspace by the index by i and denoting by J-'j 
the corresponding eigenvalue, for each state we can choose either to perform 
the Bogoliubov transformation by using the solution J^j = J^j or to leave the 
subspace unchanged by choosing J-^ = 0. At zero temperature and chemical 
potential the first choice minimizes the vacuum energy, but increasing the 
chemical potential because of Pauli blocking for an increasing number of 
states we must make the second choice. This is the mechanism for chiral 
symmetry restoration found in a four-fermion interaction model OE], and 
confirmed for gauge theories in a forthcoming paper 
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The effective mesonic action and therefore the saddle-point approximation 
can be obtained also by a variational calculation P] in which we assume as 
a test fermionic state the quasiparticle vacuum \J^). We then fix the gauge 
according to t/o(t) = 1. In the presence of such a gauge fixing we must impose 
the Gauss constraint in the Hilbert space. But since the state \J^) satisfies 
the Gauss constraint by construction, we do not need to think about it any 
longer. Under such conditions the remaining gauge fields are independent of 
time and therefore automatically satisfy periodic boundary conditions in the 
time direction. 

4.1 Background field and vacuum properties 

It might at first sight be puzzling that the form of the saddle-point solution 
does not depend on whether the theory is or is not Abelian. This point 
requires some discussion. 

We first observe that the saddle-point equations are identical to the con- 
ditions 

jf'i) = /f'2) = 0. (4.14) 

Therefore in the saddle-point approximation the fermion-antifermion mixing 
disappears in the quasiparticles action: the Bogoliubov transformations (13.21) 
in the saddle-point approximation factorize the transfer matrix into a term 
for quasiparticles and a term for antiquasiparticles [S]. Hence their effect is 
analogous to the Foldy-Wouthuysen transformation which separates positive 
from negative energy states in the Dirac Hamiltonian. With respect to this 
factorization there is no difference between Abelian and non- Abelian theories. 
The role of the condensate J-' is only to provide the background in which 
quasiparticles and antiquasiparticles propagate independently, and, as we 
will see in the next Section, it does not have a particle interpretation. For 
this reason we call J-' at the saddle point a background field. 

In order to proceed with our analysis we must distinguish two cases. In 
the first the vacuum is dominated by chromomagnetic fields with nonvanish- 
ing energy, while in the second the dominating fields are pure gauge fields. 
Here we expect a drastic difference beteen Abelian and non- Abelian theories, 
because we think that the first / second case is realized in the continuum limit 
of non- Aabelian/ Abelian gauge theories. Then in non- Abelian gauge theo- 
ries because of the nontrivial gauge-invariant vacuum, temporal link variables 
disappear from both the gauge-field and mesonic actions, but not from the 
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quasiparticle action: the QCD vacuum in the saddle-point approximation 
appears as a dual superconductor (not color superconductor) which expels 
chromoelectric fields altogether (in this dual Meissner effect the penetration 
length vanishes). Fluctuations of chromoelectric fields are subdominant, and 
as consequence we will see in Section [6] that quasiparticles are confined. In 
the Abelian case on the contrary, fiuctuations of gauge fields dominate and 
there is no confinement of quasiparticles. 

5 Symmetries and compensating fields 

In this Section we shall consider only transformations s associated with sym- 
metry groups which act in a unitary linear representation on the fermionic 
field 



and leave the action invariant. 

Since Bogoliubov transformations are unitary they preserve such symme- 
tries. But individual terms of the original action which are invariant, are 
transformed, in general, into terms which do not enjoy this property any 
longer, and symmetry conservation of the transformed total action is real- 
ized through compensations among such non invariant terms. Then there is 
the danger that approximations can disrupt such compensations resulting in 
effective symmetry breaking. 

This drawback can be avoided in many cases, requiring that the quasipar- 
ticle fields should transform in the same way as the original fermionic fields, 
because then invariant terms would obviously be transformed into invariant 
terms. This can be achieved by making the Bogoliubov transformations at 
each time slice dependent on time, and introducing, when necessary, com- 
pensating fields. 

We will restrict ourselves to symmetries which do not mix the components 
which propagate forward and backward in time, that is. 



ijj ^ ip' = sip 



(5.1) 



[s,P±] = 



(5.2) 



so that 



u' = su 




(5.3) 



14 



Then the quasiparticle operators change according to 

a' =s (^s^Rh^ (u - {s^Tsyv^^ (5.4) 
(/3')t = s (s^ kh^ {v^ + (s^J^s) u) (5.5) 

namely they are still defined by the action of s on a Bogoliubov transformed 
field where instead of J-" the modified matrix J^s is used. A simple way to 
preserve the symmetries as in the starting action is recovered if we require 
that the J-'-matrix changes under the symmetry transformation according to 

T' = sTsK (5.6) 

In order to enforce the above condition let us expand the matrices J-'t at a 
given time-slice in the basis of time-independent matrices ^{K) labeled by 
the indices K (which also include space): 

T, = J2^;iKMK) = i^,,^). (5.7) 

K 

Then (15. 6 p becomes 

(y,;,$) = (^„.$st). (5.8) 

The transformation of the basis matrices can be written as 

s <^{K) = J2 Skk'^{K') = {S ■ $)(K) (5.9) 

K' 

SO that 

iif[,^) = iift,S ■^) = iS^ -ift,^) (5.10) 

which is to say that it is necessary to require that the expansion coefficients 
transform according to 

^[{K) = Y,Sl,j,MK'). (5.11) 

K' 

The above construction also provides a physical interpretation of our formal- 
ism. Indeed we observe that, since we could perform a unitary transforma- 
tion with an arbitrary J-'-matrix and then an arbitrary expansion coefficients 
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<f{K), we can integrate over them with an arbitrary probabihty measure 
getting 

Zf = j d^i{ip) D[a\ a, /?*, /3] exp{-S'„e(J') - Sqp{a, (3; J')} . (5.12) 

Looking at the form f l3.15p of Sme we immediately reahze that time derivative 
terms are generated for the compensating fields. As a consequence, unless the 
basis matrices ^{K) are invariant, the expansion coefficients ip{K) must be- 
come dynamical bosonic fields. The basis matrices then acquire the meaning 
of structure functions of mesonic composites with quantum numbers K. K 
includes color for colored mesons, which should exist only in the deconfined 
phases. The choice of the basis matrices $(i^) (whose form must be deter- 
mined by a variational calculation) selects which mesons one will include in 
the calculation in a variational spirit. 

We can look for an approximation to this expression of the partition 
function by determining the minimum of the action with respect to 
The phases of the theory are determined by the solutions t\t oi the saddle- 
point equations. By construction if J^\j^ are matrices which minimize the 
action, then the rotated matrices J^'\ T' must also be minima of the action. 
Therefore either they coincide with the unrotated matrices or the solution 
of the minima are degenerate. And this accounts for the breaking of the 
symmetry. 

A perturbative expansion is realized by setting 

jr = 7l' + 5jrt^ jr = J + 5jr. (5.13) 

We assume the index of nilpotency of the structure functions appearing in the 
fluctuations bT\ as an asymptotic parameter and perform an expansion 
in the inverse of the nilpotency number that we call the nilpotency expansion. 
The fluctuations ^J-"^, describe in the nilpotency expansion interacting 
mesons of the form ii^v^vu^ which we will call, for easy reference, of J-'-type, 
but not of the form u^u^ -O^-O, which are not of J-'-type. A discussion about 
some mesons which are not of J-'-type is presented later in Section [71 An 
example of such an expansion can be found in Ref. 

We emphasize that from a mathematical point of view the new expression 
of the partition function is exactly equivalent to the original one. Note also 
that there is no double counting because the property of quasiparticles to 
annihilate the vacuum 

d,|J) = A|J)=0 (5.14) 
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can be interpreted as a compositeness condition: mesonic states are orthog- 
onal to quasiquark-quasiantiquark states. This constraint has the physical 
meaning of the condition Z = for bound states in the Lehmann spectral 
representation of composite operators [23l|2l] (see also |25l Vol 1, p. 461]), 
namely the condition required to introduce a bound state on the same footing 
as the constituents in a Lagrangian. 



5.1 Some examples 

We now give a few illustrative examples. The first one concerns fermion num- 
ber conservation in the nonrelativistic theory of many-body systems [3]. We 
include it because it shows clearly the need and the physics of compensat- 
ing fields and because it was the first application of our method. Moreover, 
the Bogoliubov transformation is, in this case, similar to that necessary in 
the relativistic theory of diquarks [22]; namely, it mixes the annihilation and 
creation operators of one and the same fermion (quark, electron, nucleon. 



a 



R^u-T^u^) (5.15) 

while in the relativistic case, see (13. 5p and (13. 6p . the transformation mixes 
the annihilation operator of a fermion with the creation operator of the cor- 
responding antifermion. As a consequence the operator J-" 

7 = u7u (5.16) 

carries fermion number 2 instead of zero. Under the relevant symmetry 
associated with fermion number conservation 

v! = su = e^u (5.17) 

the J-'-matrix transforms according to 

JT' = s* jr^t = e-2iejzr . (5.18) 

Since in this case the structure functions $(i^) can be taken invariant, we 
must require that 

^\K) = ^'%[K). (5.19) 

Namely, we need compensating fields ^{K^ with fermion number 2. These 
fields describe the low energy excitations. 
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We come back now to the relativistic cases of composites of fermion num- 
ber zero. The first one concerns the results obtained by the apphcation of 
our method to a four- fermion model [1], which at zero mass enjoys a discrete 
chiral symmetry generated by the parity transformation 

s = -(75®t5) (5.20) 

which commutes with the projectors P±. 

The interaction can be bi-hnearized by introducing an auxihary bosonic 

field o coupled to the fermions according to Eq. flA.4p . The relevant matrix 
in the basis is 

$ = 70 ® a (5.21) 

so that 

s$s^ = -<l> (5.22) 

and the compensating field must change sign under parity. This model is 
also interesting here because it shows that the compensating field is exactly 
the field o coupled to the fermions. 

The second example deals with the residual chiral symmetry with Kogut- 
Susskind fermions in the flavor basis. For zero fermion mass the QCD action 
is invariant under the continuous chiral transformations 

s = exp (^-^75®t5e^ (5.23) 

parametrized by the angle Q. 
The (T-field is 

a = ^^{-fQ ® 1)^ = u^{jQ ® 1) {)t + {) (7o ® 11) M (5.24) 

and the Goldstone pion which corresponds to the axial symmetry at m = 
is 

TT = z?/'^(7o75 ®h)ip = zM^(7o75 (^h)v^ + iv (7075 <^h)u. (5.25) 
We can write 

J^t = < + < (5.26) 

where the basis matrices are 

= 7o®ll, $^ = i7o75 ®t5- (5.27) 
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Under infinitesimal chiral transformations tliese basis matrices transform ac- 
cording to 

s^^s^ ^ ^^ + e^^, s^^s^ ^^^-e^^. (5.28) 

Tlierefore tlie compensating fields at, irt must transform in the inverse way. 

There are 15 more pions in four dimensions which can be constructed 
with the taste matrices t [261 [27]. Only those of J-'-type can be described 
by the fluctuations SJ-"^, . For our illustrative purposes it is sufficient to 
consider the Goldstone pion. 

The most important application concerns gauge invariance in QCD. Let 
us consider the case of a gauge transformation. In addition to (15. 3p . where s 
is replaced by the gauge transformation g{t,x), we have the transformation 
for the spatial link variables 

U',{t, x) = g{t, x) Uk{t, x) g\t, x + k) . (5.29) 

If we concentrate on colorless J-", the matrix J-'t will depend on color only 
through the configuration of spatial links U^.t 

J't = J'iUk^t) (5.30) 
where the matrices Uk^t are such that 

(f/M)x,y = ^y,x+k^'c(t,x). (5.31) 

In this case, by also introducing the matrices g and using the matrix multi- 
plication, (15. 6 p becomes 

J'iU') = J^{g Ug^)=g J^{U) g^ (5.32) 

and it is automatically satisfied. Therefore, as far as gauge invariance is 
concerned, no compensating fields are needed. 

In all the cases mentioned above the effective action respects the original 
symmetry term by term and the quasiparticle vacuum is invariant provided 
we perform the symmetry transformations on the fermions and on the com- 
pensating fields as well. Such a vacuum can be regarded as a condensate of 

the composites J-". We remark, however, that for Kogut-Susskind fermions 
these composites have a structure different from that of the chiral mesons, 

T ^a*^, + n*^^ (5.33) 
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and therefore the quasiparticle vacuum cannot be interpreted as a condensate 

of these physical particles. In this case we will refer to the field J-" as to a 
background field. But we should keep in mind that its fluctuations describe 
J-'-type mesons which include the chiral ones. 



5.2 Spontaneous chiral symmetry breaking 



We will discuss spontaneous chiral symmetry breaking in the saddle-point 
approximation. We will not derive any new results. Our purpose is only to 
formulate this problem in our formalism. 

To evaluate the order parameter, at fixed gauge configuration, we shall 
use the relation 

{i^i^)F ^-^log Zf (5.34) 



dm 



and our saddle-point approximation for the partition function, that is, 

(V'V'>^ = -^^-e, (5.35) 



r-> '-'me 1 

am 

which, we remind the reader, is justified only for the special gauge configu- 
rations we discussed earlier. By direct calculation we get 



(^^) p= tr_ i? 
and by substitution of J- 

{ipip) F = — 2 m tr 



hVTTIp. 



(70® 11)^^ + ^(70® H) 



m tr 



Hvmp. 



(5.36) 



(5.37) 



which can be expressed in terms of the eigenvalues /i„ of the Hamiltonian H 
by writing 

(V;V')F = -mJ]--L= (5.38) 

where the largest contribution in the sum comes only from the lowest eigen- 
values. 

In order to understand the meaning of this relation, let us consider the 
direct evaluation (ipip) f in the functional integral, that is (at least in the case 



2m 



—0 —hl + \pl 
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(5.39) 



where 

Po = 2 sin -— (5.40) 

is the lattice momentum in the time direction. The sum on k can be exactly 
performed (see, for example, [28l Appendix B]), and thus 



iil'il') F = —2m , coth 



^arcsmh.. 



(5.41) 



which, at zero temperature, that is in the limit of infinitely large time- 
direction Lq, is in perfect agreement with f l5.38p . 

But f l5.4ip is the starting point to derive a sort of Banks-Casher rela- 
tion [29J (see also [TTl pag. 119]) for static gauge configurations of the form 
we considered. Indeed, in the limit of vanishing lattice spacing and infinite 
volume, 

(V'V')f ~ -2m Vt- ~ / -rp{.h) (5.42) 

^ "-n Jo h 

where p(h) is the density of energy eigenvalues. Spontaneous symmetry 
breaking of chiral symmetry is recovered if in the limit of vanishing mass, 

-hm/^ r^p{h)^0. (5.43) 

This is the limit which is controlled by the density of eigenvalues near the 
origin after averaging over the gauge configurations. 

In the absence of gauge interaction (15.431) is exactly the Banks-Casher 
relation. In presence of the interaction with the gauge fields, in the Banks- 
Casher relation the spectrum of the Dirac operator is averaged in the full set 
of gauge configurations (see [30] for a recent numerical exploration). For our 
aim we are restricted instead to the energy operator, which does not contain 
time derivatives, in static gauge-configurations of the form we considered. In 
this case the number of dimensions is effectively reduced by one unit and 
therefore the relation could be more easily checked numerically. Signals that 
chiral symmetry is spontaneously broken within this restricted ensemble of 
gauge configurations would be an important check for the effectiveness of our 
approach. 
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6 Quasiparticle confinement 



In this Section we study the propagation of quasiparticles in the vacuum 
determined in the saddle-point approximation. We remind the reader that 
this assumes the dominance of chromomagnetic fields, so that the following 
developments do not apply to the Abelian case. 

At the saddle point chromoelectric fields disappear from the pure gauge- 
field and mesonic actions, so that the temporal link variables survive only in 
the action of quasiparticles. Or, in other words, the Gauss constraint still has 
to be implemented in the presence of quasiparticles. As we will see, this can 
be achieved by exactly performing the integral on temporal links variables, 
and leads to color confinement in the quasiparticle sector. 

The quasiparticle action at the saddle point, in the Uq = 1 gauge, if we 
distinguish the fields at initial and final times, reads 

^ Lo/s-l 

S qp = - [aQao + a%ai^ - PoPo - Pl^PI^) + {a^at - PtPt) 

Lo/s-l 

- J2 («:e^'^g-Wi-A+ie-^'^Q"V;) . (6.1) 
t=o 

The evaluation of the trace on the Grassmann variables, necessary at finite 
temperature, induces antiperiodic boundary conditions for the fermion fields 

a% — — ckq , aLf^ = —ao (6.2) 

f3U=-Po, Pl^ = -Po. (6.3) 

3 S 

The Gauss constraint can be implemented at a given time, say t = 
because it is conserved by the time-evolution. For this purpose we perform 
a gauge transformation at that time, 

a*Lo a% U'^ , a_Lo -> U alj, (6.4) 

s s ^ s 

Pl^^UP*^, Plo^PloUK (6.5) 

S3 3 3 

The integration on U will induce the constraint (for a discussion on the Gauss 
law in the transfer formalism of lattice gauge theories the interested reader 
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can see [HI]; a full discussion for the propagation kernel in the continuum is 
given in |3^). The fermion action becomes 



Lo/s-l 



t=l 

V s ^ S3 a / 

Lo/s-2 

- J2 («:e^^Q"'«,+i-A+ie-^'^g"V;) (6.6) 



and using the boundary conditions, 

Lo/s-l 



t=o t=o 

+ (a%_^e'^Q~'Uao-PoUU-'^Q~'P*Lo_,) ■ (6-7) 

This expression shows that U can be interpreted as the temporal link variable 
connecting time — 1 with the initial time. The effect of all other temporal 
link variables has been gauged away. 

We found convenient to perform the change of variables 

i-W-fLo , for t = 
a; = Yt, a,= { . (6.8) 

[7f_i otherwise 

(-6lo M for t = 

P: = S:, Pt=<, ' , . (6.9) 

[ot_i otherwise 

under which the quasiparticle action transforms into 

Lo/s-2 Lo/s-l 

t=0 t=0 

Lo/s-2 Lo/s-l 

+ 6lo_^Ulo6*- V 6t6:^,+ V 6te-'^Q~'6;. (6.10) 

s s 

t=0 t=0 
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The integral over U can be performed by using the result obtained in 
(see also |Ml pag- 44]) about the link integral. For SU{Nc) matrices 

j dUexp [Tr {KU^ + JU)] = exp jxr K cof (^-^^ | W{J) (6.11) 



where 



/oo 
dUexp [Tr (JU)] = J]c„(det J)" . (6.12) 
n=0 

The cof actor of any matrix A is defined by 

{cof A)ab '■= -JJI^——y^^aai...aM^-i^bbi...bN^_iAa-^bi ■ ■ ■ Aa^^_^bN^_i (6.13) 

SO that 

[A ■ (cof Af] .. = 6ij det A . (6. 14) 

In order to determine the coefficients c„, we first remark that, if d is the 
matrix with elements d/dJij, 

(det d)W{J) = J dU{det U) exp [Tr {JU)] = W{J) (6.15) 

because U G SU{Nc). But det J must satisfy the Cayley identity (see [35] 
for a complete discussion on these identities) 

(det (9) (det J)" = n{n + 1) ■■■ {n + - 1) {det J)""^ (6.16) 

so that the coefficients c„ are determined by fl6.15p to be 

1 sf (A^e - 1) sf (n) 
^" ~ ^ sf{n + N^-l) 

where sf (n) is the superfactorial of n, that is 



(6.17) 



n 



s{{n) ■=Y[k\ = Y[k''-^+\ (6.18) 

k=l k=l 

In our application we have an integral for each spatial site x with sources 



4"^'"^=-EKy*7^_,., (6-19) 

s ' ' 

^a.a.^^^a^ J^ay * _ (g ^O) 
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Then since in the present case Jxis nilpotent, with the index of nilpotency Nj 
equal to the number of quark intrinsic degrees of freedom, excluding color, 
the sum over n extends up to Nj. Then 

aet Jx — ^ , tai,...,ajve l7o,x,i • • • 10,x,i) t6i,...,6Ar^7-Lo ^ ■ ■ ■ ■ T^o _i ^ i 

(6.21) 

which is a linear combination of products of two color singlets at position x 
and times t — 0, ^ — 1, respectively. We see that, at this zero-th order of 
our perturbative expansion, at variance with J^-type mesons which already 
have a finite extension, only pointlike color singlets of quasiparticles fields 
can propagate. Indeed, at time t = 1 there are only color singlets of particles, 
and at time t = ^ only color singlets of antiparticles. Since color is conserved 
the transfer matrix cannot create colored states. 

This result can be obtained in a more concrete way by defining the transfer 
matrix for quasiparticles. To this end we first perform the change of variables 
7 e"'*^ Qj, 5 e'^^ 5 Q, and rewrite the quasiparticle action accordingly 

Lo/s-2 Lo/s-1 

" " t=o t=o 

Lo/s-2 Lo/s-1 

+ SLo_^ULoe'^Q5; - V 5te'^Q5;^,+ V 5t5; . (6.22) 

s s 

t=0 t=0 

Then we can write the quasiparticle partition function in the form 

Zqp^ J D[-fQ,-fo,So,So]{UL^-fo,SoUl^\Tqp\-fo,So) (6.23) 

where {Ul^jo,Sq 7o,(5o) are coherent states and 

%p = det (g"') exp (^f In(-e-^'^g) 7 + 5^ In(-e^'^Q)'^ (5) (6.24) 

is the quasiparticle transfer matrix. Integrating over U using the above re- 
sults we conclude that the Fock space of quasiparticles contains only pointlike 
color singlets. 

In conclusion, Zp contains the actions of baryons, antibaryons, and mesons, 
along with their interactions. The purely mesonic term with the smallest 
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number of constituents contains, in the absence of colored mesons of J-'-type, 
three quasiparticles and three antiquasiparticles. We notice, however, that 
different mesonic structures can be constructed in other ways: in terms of 
diquarks and antidiquarks, as shown in our next work [22j, or as bound states 
of J-'-mesons and quasiparticle-antiquasiparticles. 

7 Further developments 

The extraordinary results from lattice QCD push towards the attempts to 
try to recover pieces of information about baryonic interactions which are 
relevant to the phenomenology of atomic nuclei, and which cannot be ob- 
tained from phenomenology [SS], such as three-body forces and interactions 
between nucleons and strange baryons. Our method offers a way to attack 
these problems. Here we outline a derivation of an action for mesons and 
nucleons. 

In the study of the dynamics of baryons we perform the nonlinear change 
of variables in the Berezin integrals defined in [271ES] • For A*",, = 3 it reads 



onic structure functions [37] and ipt,x,i are color singlets from the triplets of 
Grassmann variables coupled to quantum numbers i. The ^t,x,2 are the new 
integration variables which are again odd elements of Grassmann algebras. 

The expansion is formulated in terms of mesonic and baryonic variables 
only, quarks being altogether eliminated. An application of this change of 
variables in a slightly different context can be found in ^37j. We remark that, 
within the limitations of validity of the approximation of that calculation, 
the resulting nucleon action contained a Wilson term as a consequence of the 
Wilson term for quarks. That effective action, therefore, did not require any 
additional care to prevent fermion doubling. 

For an illustration we evaluate the contribution quadratic in the baryonic 
variables in the present case, neglecting antibaryons and mesons of non J-"- 
type. We must then consider the expression 





t=o 



t=o 



(7.2) 
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and expand it to third order in both the 7* and 7 



^-1 

Zp ^ exp /" DYi\ 7] (1 + det J) J] 1 - ^ ( 



i:q It 



X 



n 

t=0 



(7.3) 



Now we can use the transformations (17.ip . obtaining the quadratic approxi- 
mation in the baryon variables 



Zf ^ exp {-Sme) J j D[r, ^] exp {-S, 



baryons ) 



(7.4) 



where J is the Jacobian of the transformation (17.11) . exphcitly obtained 
in [37|, and 



baryons 



{-Cij i^t+i,^^iA,^,j + ^l^^i M^,i,yj ^t,y,j) ■ (7.5) 



t=i 



The matrices appearing in the above equations are 



M 



n«2«3 Jli2j3 

V rn"^"!"!^! /'77^1\a2;'2 (T^^^xasbs 

Needless to say, odd powers of Grassmann variables always have nilpotency 
index 1, and therefore their action cannot be approximated by a nilpotency 
expansion. 

It is reasonable to assume that at low energy the important mesons are 
of J-'-type. Their interaction with baryons is 



(7.6) 



so that the mesons-nucleons action is 



d 



-M ^{K) 



d 



_dip{K) ' dip*{K) 

where S^e niust be expanded in powers of ip*, ip. 



M^*{K) 



(7.7) 



^ (7.8) 
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8 Summary and outlook 



In previous works we developed a method of bosonization of theories with 
fermions whose low energy excitations are dominated by bosonic modes. We 
were able to generate composite bosonic fields by transforming the action 
of any such theory into another exactly equivalent action. The transformed 
action can be studied in the framework of a nilpotency expansion, assuming 
as an asymptotic parameter the index of nilpotency of the composites. The 
leading approximation is given by saddle-point equations, which determine 
the properties of the vacuum. In the absence of gauge fields we solved these 
equations for both Kogut-Susskind and Wilson fermions. 

In the present work we considered the saddle-point equations for the case 
of gauge theories. This time we found one exact, gauge covariant solution 
only for Kogut-Susskind fermions. Such a solution is relevant provided the 
vacuum is dominated by chromomagnetic fields. Prom the fermionic point 
of view this vacuum appears as a condensate of composite bosons which, 
however, do not have the quantum numbers of chiral fields. We refer to the 
field of such condensed composites as to a background field. Fluctuations 
of this background field describe dynamical mesons which we call J^-type 
mesons. They describe chiral mesons, but also other mesons, including col- 
ored mesons which should not be observable at zero temperature and baryon 
density, because of a mechanism which we have not investigated. In such a 
vacuum live, in addition to mesonic J-'-type fields, fermionic quasiparticles 
with the quark quantum numbers. Thanks to the background field, quasi- 
particles do not have any direct coupling with anti-quasiparticles but are 
coupled by gauge-field interactions. 

We then explored some properties of such a vacuum. First we considered 
the spontaneous breaking of chiral symmetry, and got an expression for the 
order parameter which is accessible only in the nonperturbative regime, but 
could be evaluated in a standard Monte Carlo simulation. Second we con- 
sidered the quasiparticle action. The temporal link variables appear in this 
action in a peculiar way which allowed us to integrate them out exactly. The 
result is that only pointlike color singlets of the quasiparticle fields can prop- 
agate and therefore have a particle interpretation. They have baryonic or 
mesonic quantum numbers. Therefore, color is confined in the quasiparticle 
spectrum. This is a remarkable result in itself, and also because it allows us 
to introduce such color singlets as integration variables, and therefore as fun- 
damental fields, in the Berezin integral which defines the partition function. 
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using a formalism previously developed. Our approach, borne to introduce 
bosonic composites, has also provided a way to introduce fermionic compos- 
ites, and, therefore, the possibility of formulating QCD in terms of physical 
fields. In particular we outlined a derivation of a meson-nucleon action from 
QCD. 

Further investigation of our approach can proceed along several lines. An 
important issue is the study of the action of mesons. A calculation of this 
kind has already been done for a four-fermion model [3]. Its extension in the 
presence of gauge fields should allow us to tackle the problem of Goldstone 
fields, of the chiral anomaly and, of the utmost importance for us, the fate 
of colored mesons of J-'-type. If we could show that in the saddle-point 
approximation these mesons are also confined, the nilpotency expansion could 
be used to describe confinement and dynamics in QCD at the same time and 
on the same footing. 

Related to the above is the study of the effective action of baryons. 
This includes interactions of baryons between themselves and with mesons, 
among which those of J-'-type should dominate. All these interactions should 
also give the baryons a finite structure. They are instead pointlike in the 
saddle-point approximation, at variance with J-'-type mesons of particles- 
antiparticles, which are already extended objects. 

Another important issue concerns the theory at finite chemical poten- 
tial. We already have some results on this subject, which we will publish 
separately [22], but we have anticipated them in short form [9j. We derived 
an expression of the free energy whose numerical simulation is free of the 
sign problem. If we make the assumption of the standard theory of color 
superconductivity that at sufficiently high values of the chemical potential 
an expansion with respect to the gauge coupling constant can be justified, 
we get results compatible with the standard ones. This finding adds support 
to the physical relevance of the vacuum we have studied. 

We expect that by increasing the chemical potential and/or the temper- 
ature, an increasing number of components of the background field must 
be set equal to zero, according to a mechanism observed in the absence of 
gauge fields [SllH] , until chiral symmetry is recovered and color is deconfined. 
If this expectation is verified, the background field will result in a relevant 
parameter for the definition of QCD phases in the nilpotency expansion. 

Finally, we would like to emphasize that in comparing the present results 
with the previous results in the literature, it should be kept in mind that 
the contribution of dynamical fermions is crucial for the vacuum structure 
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in our saddle-point approximation: it changes the QCD vacuum altogether. 
More explicitly, one cannot compare with calculations in which the energy 
of the quarks is evaluated in the vacuum of the pure gauge-field theory. In 
particular, instantons are not dominant in the saddle-point approximation, 
where, instead, magnetic dipole condensation can occur. 
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A The matrices M, N of the transfer matrix 

In this Appendix we report the expressions of the matrices M, N appearing 
in the definition of the transfer matrix for the Kogut-Susskind and Wilson 
regularizations. Their common feature is that they depend only on the spatial 
link variables. 

We write the 7-matrices in terms of the Pauli matrices, adopting a con- 
vention different from that of Liischer |20j: 

^^={1 -1) ' ^^=( °a T) ' ^^=(? J) •• ^^-'^ 
Our 7-matrices have, indeed, opposite sign. 

A.l Kogut-Susskind's regularization 

Kogut-Susskind fermions in the flavor basis are defined on hypercubes whose 
sides are twice the basic lattice spacing. While in the text intrinsic quantum 
numbers and spatial coordinates were comprehensively represented by one 
index i, here we distinguish the spinorial index a = {1, . . . , 4}, the taste index 
a = {1, . . . , 4} and the flavour index i = {1, Nf}, while x = {t,xi, . . . , X3} 
is a 4- vector of even integer coordinates ranging in the intervals [0,Lt — 1] 
for the time component and [0, — 1] for each of the spatial components. 
We distinguish summations over basic lattice and hypercubes according to 
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The projection operators over fermions-antifermion states are 

P± = ^(a 1^7075 ®t5to)- (A.3) 
The matrix M = while N, neglecting an irrelevant constant, is 

3 

N = -2|(m + a) (70 ® 11) + ^(7o7j ® 1) [Pt^^f^ + ^^^^1'^] } (A-4) 
where cr is a scalar field and 



1 



are 



the lattice covariant derivative and[l] 



^There is a misprint in formula [U (A. 6)]. 
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Vj-^ = ^{^-Tt'uj) (A.6) 



P,(^) = ^(Il®Il±7,75®t5t,)- (A.7) 
The lattice hamiltonian H is related to by 

H^ = -N^N. (A.8) 
4 

Then [1] 

A= (2/J)-i + + (A.9) 
and using this expression, we derive 

n = e-"'H (yTTiP-H^ 

so that in the formal continuum limit 

o 

n ^H^ H (A.IO) 

both approach the same value. 



A. 2 Wilson's regularization 

The projection operators over fermions-antifermions are 

P± = i(]l±7o). 

The matrices M, are 

2 \2K J 

N = 2KB-^cB~^, 

where 

3 

K is the hopping parameter and 
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